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'"5 ' Abstract 

(^ . Dipion transitions of the subthreshold bottomonium levels T(nS') — * 

D I T{n'S)Tnr with n > n',n = 2,3, 4, n' = 1,2 are studied in the frame- 

work of the chiral decay Lagrangian, derived earlier. The channels 
BB,BB* + c.c,B*B* are considered in the intermediate state and 
realistic wave functions of T{nS),B and B* are used in the overlap 



l/^ I matrix elements. 

Imposing the 
ment, one obtair 
with experiment. 



f^ ■ Imposing the Adler zero requirement on the transition matrix ele- 

ment, one obtains 2d and Id dipion spectra in reasonable agreement 



o 

OO 

P ■ 1 Introduction 



K> ! Dipion transitions of heavy quarkonia first discovered in [T] , were further ex- 

H I perimentally studied in bottomonium T{nS) — > T{n' 8)1171 with n = 2, 3, n' = 

1, 2 by CLEO [21 [3l H [5], and forn = 4, n' = 1, 2 by BaBar ^ and Belle [7]. 

Recently a detailed analysis of dipion transitions between states with 
n = 1, 2, 3 was done by CLEO Collaboration [8]. 

On the theoretical side the first attempt of explanation of dipion spectra 
was done in P, [ini [HI [121 [13] using multipole gluon field expansion and 
PCAC, for a recent development of this model see [H]. 

However, in this approach the natural explanation can be given only to the 
{n,n') = (2, 1) dipion spectrum, while other types of spectra with a double 
peak need additional assumptions, such as the role of final state interaction 
and a resonance [151 [IBl [13 [IE] , exotic Tvr resonances [IHl [IHl [20l [21] , coupled 



channel effects [2^ [23] , relativistic corrections [21] , S — D mixing [23] . The 
role of constant term was studied in ^E\, for a recent development see [271125]. 
In the present paper we are using the formalism of field correlators [29] and 
chiral decay Lagrangian [30] developed for the dipion transitions in [311 [32] . 
In this formalism the dipion transition proceeds via BB, BB* etc. inter- 
mediate states and the total amplitude consists of two terms: A4 = a — b, 
where a refers to the subsequent one pion emission at each vertex of the type 
T — >■ BB, while b refers to the sequence of zero-pion and two-pion vertices, 
see Fig. 1 . The crucial for the calculation is the knowledge of the realistic wave 
functions of all participants. In [311 [32] the simplest SHO form was used, fit- 
ted to the realistic r.m.s. of a given state. In the present paper we are using 
the realistic wave functions of T{nS),B and B* mesons, calculated in [33] 
and being in good agreement with spectra and decay constants. To simplify 
calculations of the overlap matrix elements we are expanding realistic wave 
functions in series of SHO functions and check accuracy of expansion. 
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Fig. 1(a) Fig. 1 (b) 

Fig.l. Diagrams for pion emission from internal light quark loop. 
Subsequent one-pion emission (a) and two-pion emission (b). 



Another improvement over results of [311 [32] is that we also consider 
transitions (n, n') = (4, 1), (4, 2) and compare them with experiment [6l [7]. 

On more fundamental level we are considering not only spectra, but also 
absolute values of widths both for dipion and for pionless decays of the type 
T{nS) -^ BB, BB* + cc, B*B*. 

The detailed analysis, made in the paper, reveals that pionic and pionless 
decays are governed by distinct vertices, the first one is given by the chiral 
decay Lagrangian (CDL) [30], while the second one by the relativistic decay 



Lagrangian (RDL), being in some sense the relativistic generalization of the 
^Po niodel Lagrangian. 

Correspondingly we introduce two reasonable physical scales, Mbr ~ /tt ~ 
0.1 GeV and M^ ^ 2uj, where u: is the average energy of a light quark in the 
heavy-light meson [B meson), uj ^ 0.5 GeV. 

The resulting expressions are otherwise parameter-free and allow to pre- 
dict the TTTT spectrum in all transitions considered. As was already observed 
in [3n [32] , the form of the spectrum is defined by the only real parameter 
Tj, which was calculated in our approach with the account of the Adler zero 
requirement (this requirement serves as a kind of renormalization condition 
on amplitudes a,b). We shall use below the values of Adler-zero Improved 
(AZI) T] = rjAzi, obtained in [211 [32] and predict the vrvr spectrum as function 
of total TTTT mass Mtttt = q and angle 6 of tt"*" with respect to initial T di- 
rection, making also comparison to experiment. The plan of the paper is as 
follows. In section 2 general equations for dipion and pionless amplitudes are 
given, taken from [3T|, [32] , however modified as compared to [311 [32] due to 
appearance of two mass parameters M^r, M^ instead of one, Mhr in [SI]. We 
also make a correction to the Eqs. (53,55) of |3I], vahd for equal mass case, 
to make it suitable for realistic casqj. In section 3 we describe the calculation 
of overlap matrix elements and expansion of realistic wave functions. In sec- 
tion 4 and we show results of calculations for matrix elements and spectra. 
In section 5 our results are discussed in comparison with experiment. Last 
short section is devoted to conclusions and perspectives. 

2 The bottomonium decay amplitudes 

We start with the definition of Lagrangians for the light qq pair creation. As 
was derived in [30j, the CDL has the form (cf. Eq(39) in [31]) in Euclidean 
space-time 

CcDL = -l I d^X^{x)MhrU{x)4j{x) (1) 



where U{x) is the Nambu-Goldstone (NG) matrix 



/ ^ -)- ^ 7r+ K+ \ 

U = expi ij5-j-^ J , ^aK = V2 









(2) 



-'^the authors are grateful to Yu.S. Kalashnikova for pointing this fact to us. 



and /^ = 93 MeV. The Lagrangian ([T]) describes the qq pair creation with 
(or without) NG meson emission, and as will be shown below, M^r ~ fn- 

Another type of the qq pair creation occurs in the so-called time-turning 
trajectories (ttt) of a light quark at the boundary of the Wilson loop, when 
confinement ensures area law of the loop. In this case one can introduce an 
effective Lagrangian - the Relativistic decay Lagrangian (RDL) of the form 

jO-rdl = ~M d'^x'ijj{x)M^'ilj{x) (3) 

Here M^ will be shown to be of the order of 2uj , where uj is the average 
energy of light quark in B meson, u ~ 0.5 GeV (see Appendix 1 of [31] for 
exact values and discussion). Note, that M^ ^ M^r and for pionless decays 
one can neglect Ccdl as compared to Crdl- We shall use the strong decay 
formalism developed in [31j and write the amplitude Wnm for the bottomo- 
nium transition from the state n to the state m with or without additional 
NG mesons, (cf. Eq. (48) of 



-"'"(^) = 7/7^E#^^i^ (4) 

where 7 = ^, M is Mfe^ or M^^ depending on pion emission and J„„2n3 is 
the overlap matrix element of the n - th state of T and n2, n^ states of BB 
or BB* etc. For pionless decay one has 

/d Q ~ _ ~ 

yi23 7^^^n(cp + g)^„2(q)V'n3(q)- (5) 

Here c = -^^^ ~ l(fi is the energy of heavy quark in B meson, Vt ^ 
4.83 GeV) and 7/123 is the ratio of Dirac traces, 7/123 = ; ^ , where the 

Vnl^ 

projection operators Zi are defined in [31] and are of the order of 1, so that 
for T{n) -^ BB 

y.3 ^Z.^tr (,,|,,||,,^j (6) 

with S^ = {mq ± ^74 =F ip^li), Sq = mg± ^74 =F ipj^i. 

In Appendix 1 of [31] and Appendix 1 of the present work details of 
calculation of Zi are given. 



As a result one obtains for 1/123 



1/123 ~ Z 



2VL'^u 



uoVL 
qi[2il + uj) -pi 



u + n. 
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For the one-pion emission at the decay vertex one has instead 

jilmiP, k) = |yW^^„(cp - ^ + q)^n,(q)^n3(q " k). (8) 

Here I/I23 has the same origin as 7/123, but accounts for one pion emission 
at the vertex T(ra) — > {BB* + c.c.)tx, 

Finally, expanding the factor f/ in ([2]) to the second order, one obtains 
the overlap matrix element for the double pion emission vertex, 

Vii^-n^^niop - Y + q)V^n2 (q)^n3 (q - K) (10) 

Here K denotes the sum of pion momenta, K = ki + k2, and 2/123 is 

-(tttt) _ '^7ri-7r2 yi23 , . 

Having defined all overlap matrix elements in Wnm{E), Eq.(4), one can 
now express total amplitudes for processes with or without pion emission. 
The width of T{nS) due to the channel BB is given by the equation 

p PbbMbb f^c^iT ( M2 ^l nn\ 

Tn = 7c. -^ J a^pKnnanslP)! , 7a. = "^ (12) 

The dipion emission amplitude consists of two terms 

„X^^)(F) = ^ / V / ^'p jW„3(P'ki)^;;^^L(P>k2) , .. ^ 2) 



(2) 
^ / (Pp '^nn'2n^(P'kl,k2)J;;^„,^„,^(p) 



^J (27r)3E-E„.„.(p)-E(ki,k2) 

i^J (2vr)3 E-E„.„»(p) J ^ ^ 

The probability of the process T(n) -^ T{n')TXTx is obtained from tw^^T 
by standard rules 

x27r5(E(ki,k2) + E„. -K) (14) 

and the dipion decay width is 

r^-/) = j dw{{n) -^ {n'y-K) = j d^M\^ (15) 

where d^ is the phase space factor 

r2 I //\/f/\2 ^2 



1 (M2 + (M')2-g2)(M + M') [T——-^ -To r-7, , n 

^^ = 3MV^^ 4M3 ^VW^)^^\A^^4<rfgrfcos^ 

(16) 
with the notations 

M = M(T(n)), M' = M(T(n')), AM = M - M', 

g' = M^\ = {h + ^2)' = (^1 + 002? - Kl (17) 

The amplitude M. can be written accordingly to (fT3|) as 

■M=l-—j Ml 7^— -^2 (18) 

At this point we still do not impose on Ai the Adler Zero requirement. 



3 Calculation of matrix elements 

For the wave function we use the solution of the relativistic Hamiltonian, 
described in [33], (for a review see the last ref. in [29]), where the only input 
parameters are current quark masses, string tension and as- As a result one 
obtains the bottomonium spectrum with accuracy of the order of 10 MeV 
and a good agreement with experimental lepton widths, see Table 3 in [3T] . 
In what follows we call this function "the realistic wave function", meaning 
that it is among the best existing ones, but the point-by-point accuracy of it 
is not actually known. 

In this section we describe the method of calculation of 7Wi, M.2 based 
on the expansion of wave function in the full set of oscillator functions. This 
allows us to do integrals in the overlap matrix elements J, J*-^-* , J'-^-' analyt- 
ically, while the specifics of wave functions is represented by the sequence 
of numbers -coefficients in the expansion, found by the "chi squared" pro- 
cedure, namely for any radial excited state wave function found in [33] we 
write 

"'max 

^{nS,r)=Y.ct^ip,{Pr) (19) 

fc=i 

where (pk{l3r) is given in Appendix 2, and oscillator parameter (3 as well as 
coefficients c^ are obtained minimizing x^. The quality of approximations 
for different fcmax can be seen from the Fig. 2 and 3. We also compare realistic 
B meson wave functions computed in [32j with the one-term representation 
f in?]) in Fig. 4. (Note, that the B* wave function is the same as for B meson in 
the first approximation used below). Hence, keeping k = 1 for B, B* mesons, 
the simple overlap matrix element ([5]) without 1/123 is 

4°)i(p) = JtA-.Y. ci"V;.(/5i,q + cp)(^?(/32,q) = e-^/„,n(p) (20) 



fe=i 



where A = 2/5^ + /3|, and /„^ii(p) is a polynomial in powers of p^, given in 
Appendix. 

In a similar way one calculates J^^^ , J^^-* , 

_P!_k' 

4;)i(p,k)=e " ^/.,n(p)yg (21) 

_p!_k' 

41(p,ki,k2)=e " ^ «/„,n(p)yS23V.. (22) 



In (1221) we take into account that I/123 contains qi and Pi and therefore the 
result of integration, over d^q leads to the modification of /„^ii(p) = (°)/„ii(p), 
/„ii(p) -^ (^)/„ii(p). All these expressions are given in Appendix 2. 

Finally the matrix element A4 in (ITS]) can be rewritten as 

-^ = ^"p (-^^j (77J •^' - ^^ [-vi) -ir^' '''' 

and M.i,M.2 are given by the expressions 

+ (^r)'[^°^Cn'(^i)+ ^°^OM]; (24) 

M2 = zl[ (^)/::„,(o) + (i)/:„„K^i + ^2)]+ 

+ (Zr)^[ (^^Cn'(O) + ('^0(^1 + ^2)]. (25) 

Here ^''^Cf^,^, is the integral as in (IT3l) . and s = (),*, ** marks three channels: 
BB, BB* + cc] B*B* respectively; k = 0,1. Exact form for C is 

= "/„ r-pH:7) P"' 

Here v, Po,t^^\uj) and AM*^*^ are given in Appendix 3. 

Coefficients Z^ and Z2 define the relative weight of channels s = 
BB,BB* + c.c.,B*B* with account of spin and isospin, it coincides with 
the corresponding coefficients found in Table IV of ^34j- We obtain 



(Zl)^ = {Zl*f = 1, {Z;*y = 7Zl (Z*)2 = AZl Zl = -^ (27) 



PI 



Now we can compute from ^^^In,ii the decay width of the T{nS) into 
BB, BB* + c.c, B*B* etc. From ( fT2l) . taking into account, that J^«n2n3(p) = 

p2 

e-— (i)j„^„2„3(p)^, and (piPk) = ^Sikp"^, one obtains 

T{T{nS) - (5fi)W) = (^)' ^(^^^^)^e-*( «/„,.,„3(P.))^ (28) 

8 



where channel index s = {) for BB, Mj^ = ^ ^^^^^, and for s = **, decay 
channel is B*B*, M*^*^ = M%, {Z^^f = 1(Z(*)) = 4, (Z^**))^ = 7. 

4 The Adler-zero improvement of matrix el- 
ements 

As was discussed in detail in [3T] the general two-pion amplitude, consisting 
of one-pion vertices as in TWi and two-pion vertex in A^2, satisfies the Adler 
zero requirement, i.e. vanishes for kj = cjj = 0, if these amplitudes are 
imbedded in the general background not distinguishing one- and two- pion 
vertices. In particular, this implies summing up all closed channels of the 
type BB. In [31] the AZI was realized representing the amplitude A4 in Eq. 
(23) in the form 

' k? + kn ^ ^ / K2^ 

^^)a(^i,^2)-expf--. 



M = exp ( Vs2~^ 0'{^i,^2) - exp --7-52 ^(^1' ^2) (29) 



'Mbr\\, , MbrM.., 



■^1. b = ^^^^M 



2 



r j ■ '±, - n2 

\ J IT / J TY 

and requiring that a{uji = 0,^2) = b{uji = 0; U2), and the same for vanishing 

k2 = U2 = 0. 

We shall now apply the AZI procedure to our decays under consideration, 
n > n',n = 4,3,2. Before doing that, we introduce convenient new and 
universal variable x, which is in the interval [0,1] for all transitions. 

x = ^-^, f^' = {AEf - Ami (30) 

In terms of variables x, cos 6 the dipion decay probability can be written 
as 



Jo 



x(l — x) , dcos6 , , ^ , ^, ,9 



4m,| 
1^ 



-dx^-— |A<(x,cose)p. (31) 



Here Cq = ^2^3^^ =1.12-10 ^. Writing two exponents in (128!) in terms 
of x, cos 6', one has (cf.Eq. (90) of [31]) 



^ll_AE xjl-^ 

2 \ M ^ ' x + ^™' 



kf + k^ = a + 7 cos^ ^ = ^|l-^^(l-x)+ \ ^J cos^ e I (32) 

1? 
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K^ = /i^(l - x)(l - 5) = /.^(l - a;) (^1 - ^j (33) 

As it was found in [31j, the AZI amplitude (128|) vanishes for some value 
X = T], depending on the channel (n, n') and this value of rj was calculated 
in [31]. Below we shall use these values and similarly to [31] represent the 
amplitude A4 {x, cos 9) as follows 

M{x,cos6) = M{x = r],cos9) + M'{x = r],cos9){x - rf) + ... (34) 

Here prime denotes the derivative in x. It is important to note, that as 
explicit calculations show (see below), both a and b in fl28l) do not depend 
appreciably on x, and all x dependence in Ai is coming from the exponential 
factors. Using ( l32l) . ( l33l) . one arrives at the following representation for JH, 



'Mbr\ /i' 



2 



A^(x,cos0)=A^i(r7,cos^)l-^l ^e ^^^ (x-r/)/(r/), (35) 

where a = a{x = rj), 7 = 7(x = r/) and /(?]) appears due to derivatives in 
X of iD, ([33D 

AE , x2 + (l-2a;)^ 

Thus all calculations of the dipion spectra in the AZI scheme reduces to 
the calculation of A^i given in Eq. ( 12^ . at some intermediate point x = t], 
and subsequent integrations as in Eq. (!3T|) . 

This procedure refers to all {n,n') transitions of our set, except for the 
(3,2) transition. In the letter case the value of r^ = tjazi found in [31], is 
large and negative, r] = —3 and expansion (l34l) does not make sense for x 
in the physical interval [0,1]. Therefore one can instead explicitly compute 
Ml and M2 in ([23]) as given by ([21]), ([25D and insert M^ = y^u ^ 0.8 GeV 
(to be checked below by pionless decays). Computed in this way amplitudes 
M.i,M.2,V ^^^ resulting values of F^^"'" ^ = J dwT,T,{n,n') are given in the 
Table 1. 
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Table 1 

Parameters of {n,n') transitions and resulting widths rT^T^{n,n') in compari- 
son with Fexp taken from j35j . 



in,n') 


2,1 


3,1 


3,2 


4,1 


4,2 


fi, GeV 


0.483 


0.85 


0.174 


1.083 


0.479 


/i2, GeV2 


0.234 


0.721 


0.03 


1.172 


0.229 


Ml, GeV-i 


-1.56 


0.592 


-1.66 


-0.452 


0.688 


M2, GeV-i 


-0.122 


-0.116 


-0.322 


-0.0199 


-0.354 


V 





0.56 


-2.7 


0.3 


0.61 


r/(^)\keV) 


0.132 


0.195 


0.311 


0.552 


0.0071 


rexp(keV) 


6 


0.56 


0.9 


<2 


<6 



We turn now to the dipion spectra as functions of x and cos 6. For the 
AZl decay amplitudes for the processes {n, n') = (2, 1), (3, 1), (4, 1), (4, 2) one 
can use the form fl36|) . however for the 2d plots in x,cos9 one should take 
into account, that rj depends on cos 6'. Indeed, in the general AZl form of 
the dipion decays amplitude fl30l) one can see, that the first exponent on the 
r.h.s depends on cos 6* (cf. Eq. fl55]) ) hence vanishing of Ai occurs at some 
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X = r]{cos9), where 



1 (1 — ri)ri cos"^ 9 



fJicosO) =ri--fcos 9 = 7]-- ^^^ . (37) 



M" 



Correspondingly, one should replace in fl36|l r] hj fj given in fl371) . however 
this amounts to a small correction for all four transitions under investigation. 

The amplitude Ai (x, cos 9) can be expanded in a complete set of poly- 
nomials in the region < x < 1, —1 < cos 6' < 1. One can write for this 
purpose a product of orthonormal polynomials pn{x)pi{z), where p„(x) = 



\/2r7~+TP„(2x — 1), Pi{z) = J^^Pi{z), and Pkiu) is the Legendre polyno- 
mial, so that 

dw I °° 

-— = fi^Co^Jx{l -x)\M\'^ = Y. anipn{x)piicos9) (38) 

aqa cos y ^^^^ 

Qnl = dx dz—- ■^Pn{x)pi{z). (39) 

Jo J-1 dqd COS 9 

Tables of coefficients Qni for our A^(x,cos6') computed according to Eq. 
( I39l) with 7] ^ fj are given in Appendix 4. 



5 Results and discussion 

We start with the BB widths, which define our only input parameters M;,^ 
and M^. As shown in ( l30l) . for the T(45') -^ BB decay one can write 

T,siBB) = i^-^j 0.0033|4^1(p)p GeV (40) 

where JBsip) = ^xp ( — ^) hu{p), and A = 2(3f + /3| and the overlap 
integral /411 is known from the wave functions of T(45') and B meson, which 
are fitted by oscillator functions with /5i and P2 respectively. Here p = 0.33 
GeV is the B meson momentum. Approximating T(4S') wave function with 
4 and 15 oscillator functions one obtains j]j^(0.33) = —0.8434 and -3.63 
respectively, which yields T^siBB) = 2.4 MeV and 44 MeV. 

This sensitivity of r45 to the wave functions is typical for all results of 
total width, both with pion emission and without. One can exploit this fact 

12 



and invert the procedure to find the S wave function (e.g. position of zeros) 
from the values of widthqj. As it is, we consider the 15 term approximation 
for T{4S) wave function good enough (see Fig. 3 for comparison) and can 
define M^ comparing Ttheor with rcxp(T(45) -^ BB) = (20.5±2.5) MeV |35], 
which yields ^ ~ 1.46 or M^ ^ 1.72 GeV with accuracy of ±10%. However, 



the accuracy of the wave function calculation in [33] and its simulation by 
oscillator functions is not well known, and to be on the safe side, we conclude 
that Mi^ ^ 2uj as an order of magnitude estimate. We now turn to the 
total dipion widths TTjT^{nn') given in Table 1. Comparing two last lines 

in the Table for Ttheor and T^xp one can see that one gets correct order of 

r / M \ '^ 

magnitude for , "'°°;4 , so that the factor I -^ can be taken in the interval 

1 < {jtT - ^' ^-^^ 1 ^ ^ ^ 1-45- A large discrepancy for 1(21) in Table 
1 possibly is due to poor approximation of T(2S') wave function. 

In this way we are supporting our assumption, that string decay with and 
without pion emission is governed by two different scale parameters M^r and 
Mi^, which differ by one order of magnitude, M^ ^ 6-rlOMf,r. Moreover each 
of the parameters is defined by its dynamical mechanism: M^ ^ 2uj due to 
time-turning trajectories, and M^r ~ f^ due to chiral decay Lagrangian. 

Now we turn to the dipion spectra. For each transition {nn') we show 
in Figs. 5-15 three spectra:^, ^f^^g and -j-j^^, where the Id spectra are 
integrated over another variable. For comparison the existing experimental 
spectra are given in Fig. 16-18, taken from O [3, E] with theoretical curves 
from [211 [22] where rj is independent of cos^. 

One can see a general qualitative and semiquantitative agreement between 
theoretical and experimental spectra. One should stress again, that while the 
total dipion width is sensitive to the form of wavef unctions, the qualitative 
form of the spectrum is defined by the value oi rj ^ fj, which accumulates 
the information on overlap matrix elements and stabilized when the AZI 
condition is imposed. The case of the (3.2) transition has a special feature 
of low available phase space, and therefore the spectrum as function of x is 
defined mostly by the phase space factor [x(l — x)]^/^, while dependence on 
cos^ is also weak due to small /x^ = 0.03 GeV^. 

Thus all features of spectra and total widths can be understood in terms 
of the formalism, presented in [3T| |32] and in this paper above. Note, that 



^The authors are grateful to M.V.Danilov for this suggestion 
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our more refined analysis in the present paper has required two important 
changes as compared the formalism in [211 E2]- Firstly, two independent 
decay scales Mf,r and M^ are introduced and estimated here in contrast to 
the only parameter M^r in [311 [32]. Secondly, the correction was introduced 
in the P-wave vertex, viz. I/123 in Eq. (jTj), which decreases the effective value 
of this vertex. However, in the AZI form the P-wave vertices (the amplitude 
M.2 in (ESD) are derived from the Adler zero condition, as in (I25 1 ) (I551) and 
the resulting form of the total AZI amplitude is the same as in [311 [32j (with 
additional weak cos 6* dependence in fi{cos9) neglected there). 

6 Conclusions and outlook 

We have derived the amphtudes, spectra and total widths of dipion emis- 
sion in all subthreshold processes T(nS') -^ T{n' 8)1171, with n = 2,3,4 and 
n' = 1,2. We have shown that the formalism of pair creation, based on two 
effective Lagrangians. CDL and RDL with two mass parameters, M^r and 
M^ respectively, can successfully describe the experimentally found spectra. 
Moreover, it was shown that two distinct values of widths, 
Tns{BB) = 0(10 MeV) and TT,T,{nn') = 0(1 keV) can be explained by much 
different scales of M^ ^ 2uj = 0{lGeV) and M^^ ~ A = 0(100 MeV). 
Fixing in this way M^ and M^^ one obtains parameter- free one- and two- 
dimensional spectra in good qualitative and semiquantitative agreement with 
experimental data. 

The application of this formalism to the case of T(55') states and di- 
pion transitions to T{2S),T{1S) is straightforward and is now under in- 
vestigation [36j. The main difference is that the 55" state is above the 
BB, BB* + c.c, B*B* thresholds and hence the amplitudes A^i, AI2 acquire 
large imaginary parts, which strongly deform the spectra and are affected by 
final state vrvr interaction. 

Another direct application of formalism to the one-pion or rj transitions 
is clearly visible and can be a good check of the approach, planned for the 
future. 

The authors are grateful to many useful remarks and suggestions of par- 
ticipants of seminars at ITEP (Moscow) and ITP (Heidelberg). 
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Appendix 1 

Calculation of the coefficient ?/i23 in Eq. (^ 

General formalism for the calculation of the overlap matrix elements and, 
in particular, of the factor 1/123 in ([5]) is given in Appendix 1 of [31]. Here we 
present the corrected form, where all relative momenta are properly defined. 
One starts with the calculation of the trace, Eq. ([6]), which gives (cf. Eq.(52) 
of EH). 

'"^'^P.-P^^ + ^iP^.-PPoJ + ofl,)) (Al.l) 



For the (Qq)) system one can define total and relative momenta as 

Pi = Pq + Pg, qi = — ^—7^ — - (A1.2) 



Ug + Qc 



and the same for the (Qq) system. 



UJq + \Lq 



In the total cm. system Pi = p, P2 = — P and 



P? = -^1 + ^T^P ^^^-4) 

^ u + il 



Pq = -q2 -p:P; Pq = P - Pq, Pq = Pq = -P - Pq (Al.5) 

to -\- iL 



Finally one obtains from (jAl.ip . taking into account that q2 = — q^ 



1/123 



2Q^uj 



qi{2il + uj) -pi 



u + n 



(A1.6) 



The appearance of the term 0{qi) in (1A1.6P leads to the change of the 
result of integration over (i^q in Eq. ([5]) as compared to the standard inte- 
gral ^^^In,ii{p), given in (]A2.9p . Separating the factor '^, one arrives at the 
expression 

/do 

(A1.7) 
Here ak, Pi, (32, ^n,y are defined in Appendix 2, and c^ are coefficients 
of expansion of ip{nS, q) over oscillator functions, see Appendix 2 for more 
details. 



Appendix 2 

The SHO basis functions {ipk{Pr)} can be written as 

ifk{P,r)=Ck e 2 A2.1 

pr 

With the normalization condition J ipk{z)ipk'{z)d^z = 6kk' and coefficients 

_ 1 

""^ ~ (22V/2(2A;-l)!)i/2 ^^2.2) 

One can write Hermite polynomials in (]A2.ip as 

^^ = (-ir'2^;^*(-("-l).|.-^) (A2.3) 
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where 

Any radial excited state nS, n = 1,2, 3, ... can be expanded as 

oo 

V^(n^;Ar) = ^4"Vfc(/9r); (A2.5) 

fc=i 

the overlap matrix element can be written as 

Jn-xi{p) = / (0V^(^^5 q + p)^L(i^; q) (A2.6) 

where the heavy light meson wave function is 

^HLin'S; P2, r) = E €'^^^2, r) (A2.7) 



k 



fc(fci,fc2) ^^' 

X <^fci(/52,q)v^fc2(/32,q) (A2.8) 



vr 



<^fc(/3i, q + p)<^i(/52, q)v^i(/92, q) 






where 



A.^2,?+^2^ .^S^, f-^^^^ 



2/52 +/5r A„^' 



^^^/M'^' (20F//32)^ 



/?i; (22"7r3/2(2n-l)!)V2 
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Appendix 3 

CoefRcients of the overlap integrals ^^^C^I,{ijj) 

For the given intermediate state s = (),*,** thresholds Eg are 2Mb, Mb + 
Mb* , 2Mb* respectively and one can define AM*^'^^ = Eg — M, where M is 
the mass of T(nS). The resulting values of AM„^, are given below in the 
Table 2. We also assemble here some formulas for the quantities appearing 
in Eqs. ((25]), ([26]). Namely 

^-^#^' ^0 = 7^^' An = 2/?? + /?2, An' = 2/?; + /3,^ (A3.1) 

Where /3i, P[ are oscillator parameters found from fitting of realistic (nS), (n'S) 
wave functions with series of oscillator wave functions. Both /3i and f3[ de- 
pend on number k^s.x of functions kept in the series. 



2M(^) 



tW(cu) = — -^(AMW - tu), s = 0, *, * * . (A3.3) 

Po 



Other useful relations are 



kl + kl = a + -icos'e = ^ |l-^(l-x) + ^L^cos0) (A3.4) 



2 \' M ^' ' ' x + ^ 



(ki + k2)^ = /i^(l-a:)(l-^) (A3.5) 

(\ 1/2 
— — nn/ cos^ A3.6 

x+^ j ^ ' 

Table 2 

Values of mass parameter in (123]) - (126]) . 
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{n/n') 


21 


31 


32 


41 


42 


/i(GeV) 


0.483 


0.85 


0.174 


1.083 


0.479 


/i2 GeV2 


0234 


0.721 


0.03 


1.172 


0.229 


AE, GeV 


0.56 


0.895 


0.332 


1.12 


0.556 


AMnn' GeV 


-0.54 


-0.205 


-0.205 


0.02 


0.02 


AM:„, GeV 


-0.58 


-0.25 


-0.25 


-0.026 


-0.026 


AM-, GeV 


-0.625 


-0.295 


-0.295 


-0.072 


-0.072 



Appendix 4 

CoefRcients of the 2d expansion of the probabihty 



As is given in fl55]) . (15^ . the full probability 



dw 



products of polynomials 

dw 



dqd cos I 



^max. t; 



max. tin ax 



, , /, = H ainPn{x)pi{cOS( 

dqd cose „^o 



can be expanded in 



(A4.1) 



with 



Pn{x) = V2^rTlPn{2x - l),pi{z) = J'^i±lpi{z). 



(A4.2) 



The coefficients a„; calculated for ^ ^^^g in (1551) . for transitions {n,n') = 
(2, 1), (3, 1), (4, 1), (4, 2) are given below together with a„; for the (3,2) tran- 
sition, calculated as in fl23l) . Note that a;„ = for odd I. 



Table 3 The coefficients aj„ defined as in flA4.ip for the transition (2,1) 
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n 
1 





1 


2 


3 


4 





0.3890825 


0.2707011 


-0.0407149 


-0.0996852 


-0.0427755 


2 


-0.1300538 


-0.0895548 


0.0148794 


0.0335517 


0.0138329 


4 


0.0109938 


0.0072973 


-0.0016296 


-0.0029011 


-0.0010309 
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Table 4 The same as in Table 3 for the transition (3,1) 



n 
1 





1 


2 


3 


4 





0.2054341 


-0.0268913 


0.1136786 


0.0082049 


-0.0731623 


2 


-0.0739902 


0.0651164 


-0.0402644 


-0.0209861 


0.0284658 


2 


0.0145838 


-0.0245765 


0.0033265 


0.0089023 


-0.0048144 



Table 5 The same as in Table 3 for the transition (3,2) 



n 
1 





1 


2 


3 


4 





4.4761992 


0.1298668 


-1.2300728 


-0.0479114 


-0.1895183 


2 


-0.0037809 


-0.0002265 


0.0005379 


0.0001151 


0.0003569 


4 


-0.0248627 


-0.0007194 


0.0068405 


0.0002649 


0.0010495 



Table 6 The same as in Table 3 for the transition (4,1) 



n 
1 





1 


2 


3 


4 





0.6245209 


0.5257702 


0.0810291 


-0.1920402 


-0.1252025 


2 


-0.0786476 


-0.0707901 


-0.0530417 


0.0376144 


0.0335737 


4 


-0.0212058 


-0.0278274 


0.0141198 


0.0064219 


-0.0028386 
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Table 7 The same as in Table 3 for the transition (4,2) 



n 
1 





1 


2 


3 


4 





0.0165037 


-0.0083040 


0.0080081 


0.0029911 


-0.0055958 


2 


-0.0063807 


0.0055190 


-0.0024492 


-0.0020172 


0.0019898 


4 


0.0010072 


-0.0012899 


0.0000858 


0.0004939 


-0.0002193 




Fig. 2. Realistic w.f. of T(2S') from [33] (dashed line) and series of oscillator 

functions (19) with /cmax = 2. 
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Fig. 3. The same as in Fig.2, but for T(4S') and /cmax = 4. 




4 6 8 10 

r 



Fig. 4. The same as in Fig.2, but for B meson with fcmax = 1- 
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ex 

T3 




0.0 0.2 0.-1 0.6 

X 



dw{21) 



Fig.5. Dipion spectrum -^ — - as function of a; for (2,1) and (3,1) transitions. 




cos(_i>; 



Fig.6. Dipion spectrum !"^ I as function of cos 6* for the (2,1) and (3,1) transitions. 
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Fig. 7. The 2d dipion spectrum , ^ i as function of x, cos 6* for the (21) transition. 



dqd cos 




Fig.8. The same as in Fig. 7, but for the (3,1) transition. 
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M2*10 



ex 



5 



M2*6.3 



0.0 0.2 0.-1 0.6 



0.8 



Fig. 9. The same as in Fig.5, but for the (3,2) transition. Numbers on the 
curves (6,3) and (10) refer to the value of ratio -^ used in Eq.(18). 



OT 
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O 0.10 
T3 



T3 



M2*10 
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0.00 
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cos('(9) 



Fig. 10. The same as in Fig.9, but for -r-^the (3,2) transition 
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Fig.ll. The same as in Fig. 7, but for the (3,2) transition and -p^ = 10 



2.0 



T3 



T3 




0.0 0.2 0.-1 0.6 0.8 1.0 

X 



Fig. 12. The same as in Fig.5, but for the (4,1) and (4,2) transitions. 
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(4,1) 



'0.2 



O 



5 

T3 



(4,2)*10 



-1.0 -0.5 0.0 0.5 



cos(i!^) 



Fig. 13. The same as in Fig.6, but for the (4,1) and (4,2) transitions. 




Fig. 14. The same as in Fig. 7, but for the (4,1) transition. 



31 




Fig. 15. The same as in Fig. 7, but for the (4,2) transition. 
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Fig. 16. Comparison of data from [6] with parametrization (35) (broken line). 
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Fig. 17. Comparison of data from [8] with parametrization (35) (broken hne) 
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Fig. 18. Comparison of data from [6] with parametrization (35) and parame- 
ter 7] from Table 1 (dashed hne) (from |31]. 



